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Ph. The heat kernel associated with an elliptic second-order partial differential op- 

erator of Laplace type acting on smooth sections of a vector bundle over a Rieman- 
nian manifold, is studied. A general manifestly covariant method for computation 
of the coefficients of the heat kernel asymptotic expansion is developed. The tech- 
nique enables one to compute explicitly the diagonal values of the heat kernel 
coefficients, so called Hadamard-Minackshisundaram-De Witt-Seeley coefficients, 
as well as their derivatives. The elaborated technique is applicable for a manifold 
of arbitrary dimension and for a generic Riemannian metric of arbitrary signature. 
It is very algorithmic, and well suited to automated computation. The fourth heat 
kernel coefficient is computed explicitly for the first time. 

The general structure of the heat kernel coefficients is investigated in detail. 
On the one hand, the leading derivative terms in all heat kernel coefficients are 
computed. On the other hand, the generating functions in closed covariant form 
for the covariantly constant terms and some low-derivative terms in the heat kernel 
coefficients are constructed by means of purely algebraic methods. This gives, in 
particular, the whole sequence of heat kernel coefficients for an arbitrary locally 
symmetric space. 
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1 Introduction 

In this talk I am going to report on recent progress on developing some computational 
methods for the heat kernel that turned out to be very powerful for carrying out explicit 
computations Jl], 0, [| |], [7], |5], [|. We will start with the definition of the heat kernel and, 
then, will try to explain the main ideas of our approach and present the main results 
without going much into details. Consequently, I will not be as formal as probably many 
matematicians would like me to be. 

The heat kernel proved to be a very powerful tool in mathematical physics as well as 
in quantum field theory. It has been the subject of much investigation in recent years in 
mathematical as well as in physical literature. (See, for example, ]T5|, [TTj, [T3], [I], ^ [7| and 
references therein.) The study of the heat kernel is motivated, in particular, by the fact 
that it gives a general framework of covariant methods for investigating the quantum field 
theories with local gauge symmetries, such as quantum gravity and gauge theories [[UJ . 



1.1 Preliminaries 

To define the heat kernel one has to remember some preliminary facts from the differential 



geometry [11|. Let (M, g) be a smooth Riemannian manifold of dimension d with a positive 
definite Riemannian metric g. To simplify the exposition we assume additionally that it 
is compact and complete, i.e. without boundary, DM = 0. 

Let TM and T*M be the tangent and cotangent bundles of the manifold M. On 
the tangent bundle TM of a Riemannian manifold there is always a unique canonical 
connection, so called Levi-Civita connection, V™, which is torsion- free and compatible 
with the metric g. 

Let V be a smooth vector bundle over the manifold M, End ( V) be the bundle of all 
smooth endomorphisms of the vector bundle V, and C°°(M, V) and C°°(M, End (V)) be 
the spaces of all smooth sections of the vector bundles V and End (V). 

Further, we will also assume that V is a Hermitian vector bundle, i.e. there is a 
Hermitian pointwise fibre scalar product < if,ip > for any two sections of the vector 
bundle ip,if) G C°°(M, V). The dual vector bundle V* is naturally identified with V, so 
that 

< if, i/j >= ti v (f <g> if;), (1) 

where if G C°°(M, V), and if G C°°(M, V*) and tr y is the fibre trace. Using the invariant 
Riemannian volume element <ivol (x) on the manifold M we define a natural L 2 inner 
product 

{<f,ijf) = Tr L 2(<f ® ip) = Jdvo\{x)<if,4)>= J dvo\(x)tr v (if<g)ip). (2) 

M M 

The Hilbert space L 2 (M, V) is defined to be the completion of C°°(M, V) in this norm. 
Let \7 V be a connection, or covariant derivative, on the vector bundle V 

V v : C°°(M,V) -> C OD (M,T*M ®V), (3) 
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which is compatible with the Hermitian metric on the vector bundle V, i.e. 

V < ip, ip >=< V y y?, i> > + < V, > . (4) 

On the tensor product bundle T*M ® V we define the tensor product connection by 
means of the Levi-Civita connection 

v ™v = V T * M ®l + l®V y . (5) 

Similarly, we extend the connection V y with the help of the Levi-Civita connection to 
C°°{M, V^-valued tensors of all orders and denote it just by V. Usually there is no 
ambiguity and the precise meaning of the covariant derivative is always clear from the 
nature of the object it is acting on. 

The composition of two covariant derivatives is a mapping 

v t*m®v v v . c°°{M,V) -> C°°{M,T*M ®V) -> C°°(M,T*M ®T*M ®V). (6) 

Let, further, tr s denote the contraction of sections of the bundle T*M <g> T*M <g> V with 
the metric on the cotangent bundle 

tr fl = g <g> 1 : C°°(M, T*M <g> T*M ® V) -> C°°(M, V). (7) 

Then we can define a second-order differential operator, called the generalized Laplacian, 
by 

□ = tr g V T * M ® v V v (8) 

□ : C°°(M, V) -> C°°(M, T*M ® V) -> C°°(M, T*M ® T*M ® V) -> C°°(M, V). (9) 

Further, let Q be a a smooth Hermitian section of the endomorphism bundle, End(V), 
i.e. 

< </?,<2V >=< Qv,^ > ■ (io) 

Finally, we define a Laplace type differential operator F as the sum of the generalized 
Laplacian and the endomorphism Q 

F = -n+Q. (11) 
1.2 Laplace type operator in local coordinates 

The generalized Laplacian can be easily expressed in local coordinates. Let x M , (fi = 
1,2, ... ,d), be a system of local coordinates and <9 M and dx^ be the local coordinate 
frames for the tangent and the cotangent bundles. We adopt the notation that the Greek 
indices label the tensor components with respect to local coordinate frame and range from 
1 through d = dim M. Besides, a summation is always caried out over repeated indices. 
Let 

9\iv — (d/j,, d u ) be the metric on the tangent bundle, g^ v — [dx^ , dx v ) be the metric 
on the cotangent bundle, g = detg^, T^ x be the Levi-Civita connection and be the 
connection 1— form of V y . 
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Then it is not difficult to obtain for the generalized Laplacian 

□ = <rV„V„ = g~ 1/2 (d, + AJg 1 '^^ + A v ). (12) 

Therefore, a Laplace type operator is a second-order partial differential operator of the 
form 

F = -g'»d li d v -2a' i d li + q, (13) 
where a M is a End (U)-valued vector 

a» = g^A v + l -g- l/2 d u {g l l 2 g^) (14) 
and q is a section of the endomorphism bundle End (V) 

q = Q- g^ApAv - <r 1/2 d M (<? 1/2 <r A)- (15) 

Thus, a Laplace type operator is constructed from the following three pieces of geo- 
metric data 

• a metric g on M, which determines the second-order part; 

• a connection 1-form A on the vector bundle V, which determines the first-order 
part; 

• an endomorphism Q of the vector bundle V, which determines the zeroth order part. 

It is worth noting that every second-order differential operator with a scalar leading 
symbol given by the metric tensor is of Laplace type and can be put in this form by 
choosing the appropriate connection V v and the endomorphism Q. 

1.3 Self-adjoint operators 

Using the L 2 inner product we define the adjoint F* of a differential operator F by 

(F*<p,il>) = (<p,Fil>). (16) 

It is not difficult to prove that if the connection V is compatible with the Hermitian 
metric on the vector bundle V and the boundary of the manifold M is empty, then 
the generalized Laplacian □, and, obviously, any Laplace type operator F, is an elliptic 
symmetric differential operator 

(□^V) = (y,n^), {Fip,il>) = {ip,Fil)), (17) 

with a positive principal symbol. Moreover, the operator F is essentially self-adjoint, i.e. 
there is a unique self-adjoint extension F of the operator F. We will not be very careful 
about distinguishing between the operator F and its closure F, and will simply say that 
the operator F is elliptic and self-adjoint. 
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Spectral theorem. There is a well known theorem about the spectrum of any elliptic 
self-adjoint differential operator F acting on smooth sections of a vector bundle V over 
a compact manifold M, F : C°°(M, V) — > C°°(M,V), with a positive definite principal 



symbol |15j. Namely 



• the operator F has a discrete real spectrum, A n , (n = 1,2,...,), bounded from 
below 

A n > -C, (18) 

with some real sonstant C, 

• all eigenspaces of the operator F are finite-dimensional and the eigenvectors, ip n , of 
the operator F, 

Fip n = X n Lp n , (19) 

are smooth sections of the vector bundle V, which form a complete orthonormal 
basis in L 2 (M, V). 

{Vm <Pm) = 5 mn . (20) 

In the following we wil assume that the endomorphism Q is bounded from below 
by a sufficiently large constant, so that the Laplace type operator F is strictly positive. 
This does not influence all the conclusions but simplifies significantly the technical details 
needed to treat the negative and zero modes of the operator F . This can be always 
done as long as we study only asymptotic properties of the spectrum for large eigenvalues 
but not the structure and the dimension of the null space and related cohomolgical and 
topological questions. 



1.4 Heat kernel 

Thus all eigenfunctions of the Laplace type operator F are smooth sections of the vector 
bundle V and, if the manifold M is compact, F has a unique self-adjoint extension, which 
we denote by the same symbol F. 

Then the operator U(t) = exp(-tF) for t > is well defined as a bounded operator on 
the Hilbert space of square integrable sections of the vector bundle V. These operators 
form a one-parameter semi-group. The kernel U(t\x, x') of this operator is defined by 

U{t\x,x) = exp(-tF)8(x,x') = £ e -' A "^ n (x) <g> <p n (x'), (21) 

n 

where 8(x,x') is the covariant Dirac distribution along the diagonal of M x M, and each 
eigenvalue is counted with multiplicities. It can be regarded as an endomorphism from 
the fiber of V over x' to the fiber of V over x. 

The kernel U(t\x,x') of the operator exp(— tF) satisfies the heat equation 

(d t + F)U(t\x,x') = (22) 



with the initial condition 

U(0 + \x,x') = 5(x,x'). 



(23) 
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That is why, it is called the heat kernel. It can be proved that there is a unique smooth 
solution, called the fundamental solution, of the heat equation satisfying that initial con- 
dition. Thus, the heat kernel is the fundamental solution of the heat equation. For t > 
the heat kernel is a smooth section of the external tensor product of the vector bundles 
VmV* over the tensor product manifold M x M: U(t\x,x') G C°°(]R + x M x M, VmV*). 

It is not difficult to prove that for a positive elliptic operator F the inverse operator, 
called also Green operator, G = F _1 , is a bounded operator with the kernel given by 

oo 

G(x,x') = J dtU(t\x,x'). (24) 
o 

1.5 Trace of the heat kernel and the spectral functions 

As we already said above, for alH > the heat semi-group U(t) = exp(-tF) of a Laplace 
type operator F on a compact manifold M is a bounded operator on the Hilbert space 
L 2 (M, V) and is trace-class, with a well defined trace given by the formula 

Tr L2 exp(-tF) = f dvol (x)tr v U(t\x, x) = ^e _ * A ". (25) 

M 

The trace of the heat kernel is obviously a spectral invariant of the operator F . It 
determines all other spectral functions by integral transforms. 

1. The distribution function, N{\), defined as the number of the eigenvectors with the 
eigenvalues less than A, is given by 

c+ioo 

N{\) = #{v? n | A n < A} = — J — e* A Tr L 2 exp(— tF), (26) 

c—ioo 

where c is a positive constant. 

2. The density function, p(A), is defined by derivative of the distribution function and 
is obviously 

c+ioo 

P(A) = d\ N{X) = 2ri I dtetXTY L*exp(-tF). (27) 

c—ioo 

3. The zeta-f unction, ((s), defined as the trace of the complex power of the operator 
F, is given by 

1 °° 

C(s) = Tr L 2F~ S = J d\p(\)\~ s = — J dtf- 1 Tr L 2 exp(-tF), (28) 



where s is a complex variable with Res > d/2. 
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4. Finally, the trace of the powers of the resolvent operator is also expressed in terms 
of the heat kernel 

oo 

R k (z) = Ti L 2(F-z)- k = J d\p(X)(X-z)- k = —-— ] dtt fe -V 2 Tr L2 exp(-tF), 

(29) 

where z is a complex parameter with Re z < and k is a positive integer satisfying 
the condition k > d/2. 

These spectral functions are very useful tools in studying the spectrum of the operator 
F . The zeta function enables one to define, in particular, the regularized determinant of 
the operator F, 

C'(0) = -logDetF, (30) 

which determines the one-loop effective action in quantum field theory All these functions 
are, in principle, equivalent to each other. However, the heat kernel is a smooth function 
whereas the distribution and especially the density function are extremely singular. That 
is why the heat kernel seems to be more convenient for practical purposes. 



2 Asymptotic expansion of the heat kernel 

In the following we are going to study the heat kernel only locally, i.e. in the neigh- 
bourhood of the diagonal of M x M, when the points x and x' are close to each other. 
The exposition will follow mainly our papers |], We will keep a point x' of the mani- 
fold fixed and consider a small geodesic ball, i.e. a small neighbourhood of the point x'\ 
B x i = {iG M\r(x, x') < e}, r(x, x') being the geodesic distance between the points x and 
x' . We will take the radius of the ball sufficiently small, so that each point x of the ball 
of this neighbourhood can be connected by a unique geodesic with the point x' . This can 
be always done if the size of the ball is smaller than the injectivity radius of the manifold 
at x', e < r in j(x'). 

Let a(x, x') be the geodetic interval, also caled world function, defined as one half the 
square of the length of the geodesic connecting the points x and x' 

a(x, x') = -r 2 (x, x'). (31) 

The first derivatives of this function with respect to x and x' define tangent vector fields 
to the geodesic at the points x and x' 

u » = g» u V v <7, u»' = £ mV V>. (32) 

and the determinant of the mixed second derivatives defines a so called Van Vleck-Morette 
determinant 

A(x, x') = g-V 2 (x)det (-V M V>(x, x'))<r 1/2 (x')- (33) 

Let, finally, V(x,x') denote the paralell transport operator along the geodesic from the 
point x' to the point x. It is a section of the external tensor product of the vector bundle 
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VSV* over M x M, or, in other words, it is an endomorphism from the fiber of V over 
x' to the fiber of V over x. 

Near the diagonal ofMxM all these two-point functions are smooth single-valued 
functions of the coordinates of the points x and x' . To simplify the consideration, we will 
assume that these functions are analytic. 

Then, the function 

U (t\x, x') = (47it)- d/2 A(x, x') exp {-—ofa a;')) V(x, x') (34) 
satisfies the initial condition 

U (0 + \x,x r ) = 5(x,x'). (35) 

Moreover, locally it satisfies also the heat equation in the free case, when the Rieman- 
nian curvature of the manifold Riem, the curvature of the bundle connection 1Z and the 
endomorphism Q vanish: Riem = 1Z — Q — 0. Therefore, U (t\x : x') is the exact heat 
kernel for a pure generalized Laplacian in flat Euclidean space with a flat trivial bundle 
connection and without the endomorphism Q. 

This function gives a good framework for the approximate solution in the general case. 
Namely, by factorizing out this free factor we get an ansatz 

U(t\x, x') = (47rt)- d/2 A(:r, x') exp [-^(x, x'^j V(x, x')n(t\x, x'). (36) 

The function Q(t\x,x'), called the transport function, is a section of the endomorphism 
vector bundle End (V) over the point x'. Using the definition of the functions a(x, x'), 
A(x,x') and V(x,x') it is not difficult to find that the transport function satisfies a 
transport equation 

(d t + ±D + L) n(t) = 0, (37) 

where D is the radial vector field, i.e. operator of differentiation along the geodesic, 
defined by 

D = V u = «"V„, (38) 
and L is a second-order differential operator defined by 

L = V^A-^FA^V. (39) 

The initial condition for the transpot function is obviously 

n(t\x,x') = I, (40) 

where / is the identity endomorphism of the vector bundle V over x' . 

One can prove that when the operator F is positive the function Q(t) satisfies the 
following asymptotic conditions 

lim t a d?tt(t) = for any a > 0, N> 0. (41) 

t— >oo,0 
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In other words, as t — » oo the function Q(t) and all its derivatives decreases faster than 
any power of t, actually it decreases exponentialy, and as t —>■ the product of Q(t) with 
any positive power of t vanishes. 

Now, let us consider a slightly modified version of the Mellin transform of the function 
Q(t) introduced in |2| 



— / dtr^nti). (42) 

-q) Jo 



b < = r ( -,) 

This integral converges for Reg < 0. By integrating by parts iV times and using the 
asymptotic conditions ( [IT] ) we get also 

1 POO 

\= v( - m / dtr«-i+ N (-d t ) N m- (43) 
r(— q + N) Jo 

This integral converges for Reg < N — 1. 

Using this representtion one can prove that |2| 

• the function b q is analytic everywhere, i.e. it is an entire function, 

• the values of the function b q at the integer positive points are given by 

b k = (-d t ) k Q(t)\ t=o , (44) 

• b q satisfies an asymptotic condition 

lim T(-g + N)b q = 0, for any N > 0. (45) 

\q\— >oo, Reg<iV 

By inverting the Mellin transform we obtain a new ansatz for the transport function 
and, hence, for the heat kernel 

c+ioo 

fi(t) = — / dqt«T{-q)b q (46) 
2m J 

c—ioo 

where c < is a negative constant. 

Substituting this ansatz into the transport equation we get a functional eqution for 
the function b q 

[l+ l -D^b q = Lb q ^. (47) 
The initial condition for the transport function is translated into 

bo = I- (48) 

Thus, we have reduced the problem of solving the heat equation to the following 
problem: one has to find an entire function b q (x, x') that satisfies the functional equation 
([17]) with the initial condition fl48|) and the asymptotic condition (fHf). 
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The function b q turns out to be extremely useful in computing the heat kernel, the 
resolvent kernel, the zeta-function and the determinant of the operator F. It contains the 
same the information about the manifold as the heat kernel. In some cases the function 
b q can be constructed just by analytical continuation from the integer positive values bk 

Now we are going to do the usual trick, namely, to move the contour of integration 
over q to the right. Due to the presence of the gamma function T(—q) the integrand has 
simple poles at the non-negative integer points q = 0, 1, 2 . . ., which contribute to the 
integral while moving the contour. So, we get 



fi W= E ^fh + R N (t), (49) 



N-l ( _+\k 

£ 

fc=0 



where 



cjv+ioo 



R N (t) = ^- I dqt q T{-q)b q (50) 

with Cat is a constant satisfying the condition N — I < cn < N. Asi^O the rest term 
i?jv(^) behaves like 0(t N ), so we obtain an asymptotic expansion as t — > 

n(t\x,x')~j2^rb k (x,x'). (si) 

Using our ansatz (|36|) we find immediately the trace of the heat kernel 

c+ioo 

Tr i2 exp(-tF) = (4vrt)- d / 2 -l- / dqt* T(-q)B q , (52) 

2m J 

c— ioo 

where 

B q = Ti L 2b q = J dvo\ (x)tv v b q (x, x). (53) 

M 

The trace of the heat kernel has an analogous asymptotic expansion as t — > 

Tr L2 exp(-tF) ~ (4vrt)^/ 2 £ ^B k . (54) 

fc>0 fc - 

This is the famous Minackshisundaram-Pleijel asymptotic expansion. The physicists 
call it the Schwinger-De Witt expansion |L0] . Its coefficients Bk are also called sometimes 



Hadamard-Minackshisundaram-De Witt-Seeley (HMDS) coefficients. This expansion is 
of great importance in differential geometry, spectral geometry, quantum field theory and 
other areas of mathematical physics, such as theory of Huygence' principle, heat kernel 
proofs of the index theorems, Korteveg-De Vries hierarchy, Brownian motion etc. (see, 
forexample, |nj). 



For integer q — k — 1, 2, . . . the functional equation ( 47|) becomes a recursion system 



that, together with the initial condition (|48f) , determines all the HMDS-coefficients bk- 
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One should stress, however, that this series does not converge, in general. In that 
sense our ansatz (fj9|) in form of a Mellin transform of an entire function is much better 
since it is exact and gives an explicit formula for the rest term. 

Let us apply our ansatz for computation of the complex power of the operator F 
defined by 

oo 

qp = f -p = __ j dtt P-i Tj(ty (55) 
Using our ansatz for the heat kernel we obtain 



o 



c—ioo 

where c < —Rep + d/2. 

Outside the digonal, i.e. for d / 0, this integral converges for any p and defines an 
entire function of p. The integrand in this formula is a meromorphic function of p with 
some simple and maybe some double poles. If we move the contour of integrtion to the 
right, we get contributions from the simple poles in form of powers of a and a logarithmic 
part due to the double poles (if any). This gives the complete structure of diagonal 
singularities of the complex power of the opertor F, G p (x, x'). Thus the function b q turns 
out to be very useful to study the diagonal singularities. 

In particular case p = 1 we recover in this way the singularity structure of the Green 
function 

G = (Any^A^V ($ + * log |) + G reg , (57) 



where 



(_-l\k / r > ,d/2-l-k 

*= E Hr r(rf/2 - 1 - k ) (-) b * ( 58 ) 



0, for odd d 

* = I (-l)d/2 (59) 

—. - , . bd/2-i for even d 
T(d/2) 1 

G Ieg = {^Y^^V— J dqT(-q)T(-q-l + d/2)(^j b q , (60) 

a—ioo 

where [d/2] — 1 < a < [d/2] — 1/2. We see that due to the absence of the double poles 
in the integrand there is no logarithmic singularity in odd dimensions. Thus, all the 
singularities of the Green function are determined by the HMDS-coefficients bk and the 
regular part is determined by the function b q as well. 

Now, let us consider the diagional limit of G p . By taking the limit o — > we obtain a 
very simple formula in terms of the function b q 

G p (x,x) = (4n)- d / 2 T{p ~^ /2) b d/ ^ p (x,x). (61) 
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This gives automatically the zeta-function of the operator F 



C(P) = ^r dl2 T \^ 2) B d/ ^ P . (62) 



Herefrom we see that both G p (x, x) and ((p) are meromorphic functions with simple 
poles at the points p = [d/2] + 1/2 - k, (k = 0, 1, 2, . . .) and p = 1, 2, ... , [d/2]. In 
particular, the zeta-function is analytic at the origin. Its value at the origin is given by 

] for odd d 

C(0) = { (-IY' 2 (63) 

This gives the regularized number of all modes of the operator F since formally 

C(0) = Tr L2 l = £l. (64) 

n 

Moreover, the derivative of the zeta-function at the origin is also well defined. As 
we already mentioned above it determines the regularized determinant of the operator F 
since formally 

log Det F = Tr L2 log F = £ log A n = -('(0). (65) 



Thus we obtain for the determinant 



n (-l)(d+l)/2 

logDetF = ~W /2 Y(d/ 2 + l) B * /2 ^ ° M * (66) 



and 



logDetF = (47r)- rf / 2 r( l rf/2 j +i) {B' d/2 - [*(d/2 + 1) + C)B d/2 } for even d. (67) 



Here ^(z) = (d/dz)logT(z) is the psi-function, C = — ^(1) is the Euler constant, and 

(6? 



= 



q=d/2 



3 Non-recursive solution of the recursion system 

The main problem we are solving is to compute the HMDS-coefficients, not only the 
integrated ones B k = J M dvol (x)tr vbk(x, which are determined by the diagonal values 
of bk(x,x), but rather the off-diagonal coefficients bk(x,x'). They are determined by a 
recursion system which is obtained simply by restricting the complex variable q in the 
eq. (H7I) to positive integer values q = 1,2,.... This problem was solved in f|Tj, [2j where 
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a systhematic technique for calculation of bk was developed. The formal solution of this 
recursion system is 

So, the problem is to give a precise practical meaning to this formal operator solution. 
To do this one has, first of all, to define the inverse operator (1 + D/k)^ 1 . This can be 
done by constructing the complete set of eigenvectors of the operator D. However, first 
we introduce some auxiliary notions from the theory of symmetric tensors. 



3.1 Algebra of symmetric tensors 

Let u a and e a be the basises in the cotangent T*M and tangent TM bundles, S n (M) 
be the bundle of symmetric contravariant tensors of rank n, S n (M) be the bundle of 
symmetric n-forms and S^(M) = S m (M) <g> S n (M) be the bundle of symmetric tensors 
of type (m, n) with the basis 

sll'.'.K = u(ai ® • • • ® ^ am) ® e(&! <E> • • • ® e K) . (70) 

where the parenthesis mean the symmetrization over all indices included. 
In the space S 1 ™ there is a natural unity symmetric tensor 

hn) = s{:: n n , (7i) 

which is an identical endomorphism of the vector bundles S n and S n . 
We define the following binary operations on symmetric tensors: 

a) the exterior symmetric tensor product V 

V : S™ x Sj — > (72) 

by 

A \ I T3 — / 

(ai.. .(Ira a m+l ••■O'm+j 

b) and an inner product * 

* : S:xS;- Sj,, (74) 

by: 

^^ = ^::r^:±C:.r- ( 75 ) 

Further, we define also an exterior symmetric covariant derivative sV on symmetric 
tensors 

sV : S% - S™ +1 (76) 

by 

s vA = v (ai <:t +l) < 1 :t +1 - (77) 

Everything said above remains true if we consider End (y)-valued symmetric tensors, 
i.e. sections of the vector bundle <g> End (V), for some vector bundle V over M. The 
product operations include then the usual endomorphism (matrix) inner product as well. 



y q _ j^p\...b n j^b n+1 ...b n+i ) ai...a m+ j (73) 

(a-i.. • Qirn. dm. -1-1 ■ ■ •CLrti.A-'i 
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3.2 Covariant Taylor basis 

Let us consider the space C{B x i) = {\f >= f(x, x')\ x e B x i} of smooth analytic two- 
point functions in a small neighbourhod B' x of the diagonal x = x' . Here we denote the 
elements of this space by \ f >. Let us define a special set of such functions \n >G C(B x i) 
labeled by a natural number n G N by 

|0 > = 1 

\n> = ^-y n u\ (ra = 1,2,...), (78) 
n! 

where u' is the tangent vector field to the geodesic conecting the points x and x' at the 
point x' given by the first derivative of the geodetic interval a 

v! = (g' ab V b o-)e> a , (79) 

where prime ' denotes the objects and operations at the point x'. The functions \n > are 
two-point geometric objects, which are scalars at the point x and symmetric contravariant 
tensors at the point x', more precisely, they are sections of the vector bundle S n over the 
point x'. 

Let us define also the dual space of linear functionals 

C*(B X ,) = {< f\: C(B X ,) - C}, (80) 

with the basis < n\ dual to the basis \n >. The values of the dual basis functionals on 
the two-point functions are sections of the vector bundle of symmetric forms S n defined 
to be the diagonal values of the symmetric exterior covariant derivative sV 

<n\f>=[(sV) n f\, (81) 

where the square brackets mean restriction to the diagonal x — x' . 
The basis < n\ is dual to \m> in the sense that 

< n \m >= 8 mn I(n). (82) 

Using this notation the covariant Taylor series for an analytic function \f > can be 
written in the form 

\f>=J2\n>*<n\f>. (83) 

n>0 

Now it is almost obvious that our set of functions \n > forms a complete basis in 
C(B x i) due to the fact that there does not exist any nontrivial analytic function which is 
'orthogonal' to all of the eigenfunctions \n >. In other words, an analytic function that 
is equal to zero together with all symmetrized derivatives at the point x — x' is, in fact, 
identically equal to zero in B x i. 

It is easy to show that these functions satisfy the equation 



(84) 
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and, hence, are the eigenfunctions of the operator D with positive integer eigenvalues. 
Note, however, that the space of analytical functions C(B X >) is not a Hilbert space with a 
scalar product < f\g > defined above since there are a lot of analytic functions for which 
the norm < f\f > diverges. If we restrict ourselves to polynomial functions of some order 
then this problem does not appear. Thus the space of polynomials is a Hilbert space with 
the inner product defined above. 



3.3 Covariant Taylor series for HMDS-coefficients bt 

The complete set of eigenfunctions |n > can be employed to present an arbitrary linear 
differential operator L in the form 

L = V] \m > -k < m\L\n > -k < n\, (85) 

m,n>0 

where < m\L\n > are the 'matrix elements' of the operator L that are just End (U)-valued 
symmetric tensors, i.e. sections of the vector bundle S^(M) ® End (V). We will not study 
the question of convergency of the expansion (|85D . It can be regarded just as a formal 
series. When acting on an analytic function, this series is nothing but the Taylor series 
and converges in a sufficiently small region B x i. 

Now it should be clear that the inverse operator (1 + \D)~ l can be defined by 

l + T D) =Y |n>*<ra|. (86) 

Using this representation together with the analogous one for the operator L, (^J) , we 
obtain a covariant Taylor series for the coefficients 

h = J2\ n> * < n \ b k > ( 87 ) 

n>0 

with the covariant Taylor coefficients < n\bk > given by [T|, 

k k-1 1 



< n\b k >= 



n 1 ,...,n k _ 1 >o k + n *;-l + n fe _i 1 + ni 
x < n\L\nk-i > * < rife_i|X|n^_2 >*•••*< ni|L|0 >, (88) 

where < m\L\n > are the matrix elements of the operator L (|5U|). 

It is not difficult to show that for a differential operator L of second order, the matrix 
elements < m\L\n > do not vanish only for n < m + 2. Therefore, the sum (|88f) always 
contains only a finite number of terms, i.e., the summation over rij is limited from above 

m > 0, Ui<n i+1 + 2, (i = 1, . . . , k - 1; n k = n). (89) 
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3.4 Matrix elements < m\L\n > 

Thus we reduced the problem of computation of the HMDS-coefficients bk to the com- 
putation of the matrix elements of the operator L. The matrix elements < m\L\n > are 
symmetric tensors of the type (m,n), i.e. sections of the vector bundle S^(M). 

The matrix elements < n\L\m > of a Laplace type operator have been computed in 
our papers [||, [|. They have the following general form 

< m\L\m + 2 >= g' 1 V 7 (m) (90) 

< m\L\m + 1 >= (91) 

(7Tl\ ( Tfl \ 

J /(„) V Z( m _ n ) - ( -J hn-l) V F( m _ n+ i) 

(Tfl \ 
I(n-2) V X( 
m—n+2): 
(92) 

where g~ x is the metric on the cotangent bundle, Z< n \ is a section of the vector bundle 
S n {M) eg) End (\^), Y( n ) is a section of the vector bundle S^(M) <S> End (V) and is a 
section of the vector bundle S^(M). Here it is meant also that the binomial coefficient 
m is equal to zero if < or n < k. 

We will not present here explicit formulas for the objects Z(„), Y( n ), an d X^, (they 
have been computed for arbitrary n in the our papers 0,0), but note that all these 
quantities are expressed polynomially in terms of three sorts of geometric data: 

• symmetric tensors of type (2,n), i.e. sections of the vector bundle S%(M) 

K (n) = (sV) n ~ 2 Riem, (93) 
where Riem is the symmetrized Riemann tensor 

Riem = R% d \ )S %, (94) 

• sections of the vector bundle End (V) £g> S^(M) 

K {n) = (sV^n, (95) 
where TZ is the curvature of the connection on the vector bundle V in the form 

n = K\sl (96) 

• End (V^)-valued symmetric forms, i.e. sections of the vector bundle End (V)(g>S n (M), 
constructed from the symmetrized covariant derivatives of the endomorphism Q of 
the vector bundle V 

Q(n) = (sV) n Q. (97) 
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From the dimensional arguments it is obvious that the matrix elements < n\L\n > are 
expressed in terms of the Riemann curvature tensor, Riem, the bundle curvature, TZ, and 
the endomorphism Q; the matrix elements < n + l|L|n > — in terms of the quantities 
VRiem, VTZ and VQ; the elements < n + 2\L\n > — in terms of the quantities of the 
form VVRiem, Riem ■ Riem, etc. 

3.5 Diagramatic technique 

In the computation of the HMDS-coefficients by means of the matrix algorithm a "dia- 
grammatic" technique, i.e., a graphic method for enumerating the different terms of the 
sum (|88|) , turns out to be very convenient and pictorial [ij. 

The matrix elements < m\L\n > are presented by some blocks with m lines coming 
in from the left and n lines going out to the right (Fig. 1), 




Fig. 1 

and the product of the matrix elements < m\L\k > * < k\L\n > - -by two blocks 
connected by k intermediate lines (Fig. 2), 




Fig. 2 

that represents the contractions of the corresponding tensor indices (the inner product). 

To obtain the coefficient < n\bk > one should draw, first, all possible diagrams which 
have n lines incoming from the left and which are constructed from k blocks connected 
in all possible ways by any number of intermediate lines. When doing this, one should 
keep in mind that the number of the lines, going out of any block, cannot be greater 
than the number of the lines, coming in, by more than two and by exactly one. Then 
one should sum up all diagrams with the weight determined for each diagram by the 
number of intermediate lines from the analytical formula fl5B|). Drawing of such diagrams 
is of no difficulties. This helps to keep under control the whole variety of different terms. 
Therefore, the main problem is reduced to the computation of some standard blocks, 
which can be computed once and for all. 

For example, the diagrams for the diagonal values of the HMDS-coefficients [bk] =< 0\bk > 
have the form, 



[61] = o 



(98) 
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N =00+3 XX) 



(99) 



fed 
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(100) 



3.6 Remarks 

Let us make some remarks about the elaborated technique. 

• this technique is applicable for a generic Riemannian manifold M and for a generic 
vector bundle V of arbitrary dimensions, 

• this technique is manifestly covariant, which is an inestimable advantage in quantum 
field theory, especially in quantum gravity and gauge theories, 

• since it is purely local, it is also valid for manifolds with boundary and noncom- 
pact manifolds, provided one considers the local HMDS-coeflicients bk(x,x') in a 
small neighbourhood B of the diagonal ofMxM that does not intersect with the 
boundary, B n dM = 0, 

• moreover, this technique works also in the case of pseudo-Riemannian manifolds 
and hyperbolic differential operators, 

• this method is direct — it does not need to use any sophisticated functorial properties 
of the integrated coefficients B k , 

• it gives not only the diagonal values of the HMDS-coefficients [b k ] but also the 
diagonal values of all their derivatives; thus it gives immediately the asymptotics of 
the trace of derivatives of the heat kernel 



where P is a differential operator, 

• due to use of symmetric forms and symmetric covariant derivatives the famous 
'combinatorial explosion 1 in the complexity of the HMDS-coefficients is avoided, 

• the developed technique is very algorithmic and well suited to automated compu- 
tation — there are a number of usual algebraic operations on symmetric tensors 
that seems to be easily programmed, the needed input, i.e. the matrix elements 
< n\L\m >, is computed in advance analytically and is already known, 



Tr L2 Pexp(-tF), 



(101) 
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the developed method is very powerful] it enabled us to compute for the first time 
the diagonal value of the fourth HMDS-coefficient [64] [l], [|. (The third coefficient 



[63] was computed first by Gilkey [fj|. 



last, this technique enables one not only to carry out explicit computations, but also 
to analyse the general structure of the HMDS-coefficients for all orders k. 



4 Covariant approximation schemes for the heat ker- 
nel 

4.1 General structure of HMDS-coefficients 

Now we are going to investigate the general structure of the HMDS-coefficients. We will 
follow mainly our papers [p], @, 0, f|, ||, || (see also our review papers |7|, [| 3)- Our 
analysis will be again purely local. Since locally one can always expand the metric, the 
connection and the endomorphism Q in the covariant Taylor series, they are completely 
characterized by their Taylor coefficients, i.e. the covariant derivatives of the curvatures, 
more precisely by the objects K^, TZi n ) and Q^ n ) introduced above. We introduce the 
following notation for all of them 

= {K (n+2) ,K {n+1) , Q ( „)}, in = 0, 1, 2, . . .), (102) 

and call these objects covariant jets, n will be called the order of a jet Further we 

introduce an infinite set of covariant jets of all orders 

.7 = {&(»); (n = 0,1,2,...)}. (103) 

The first two HMDS-coefficients have a very well known form [ 14|, § 



B = / dvo\{x)tr v I. (104) 



M 

B x = j dvol (x)tr v (q - ^Rj , (105) 

M 

where R is the scalar curvature. 

As far as the higher order coefficients B^, (k > 2), are concerned they are integrals 



of local invariants which are polynomial in the jets ||15| . One can classify all the terms 
in them according to the number of the jets and their order. The terms linear in the 
jets in higher order coefficients B^, (k > 2), are given by integrals of total derivatives, 
symbolically f M dvol (x) tr v □ fc ~ 1 3ft. They are calculated explicitly in |L]. Since the total 
derivative do not contribute to an integral over a complete compact manifold, it is clear 
that the linear terms vanish. Thus Bk, (k = 2, 3, . . .), begin with the terms quadratic in 
the jets. These terms contain the jets of highest order (or the leading derivatives of the 
curvatures) and can be shown to be of the form f M dvol (x)tr y 3? □ k ~ 2 3l. Then it follows a 
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class of terms cubic in the jets etc.. The last class of terms does not contain any covariant 
derivatives at all but only the powers of the curvatures. In other words, the higher order 
HMDS-coefficients have a general structure, which can be presented symbolically in the 
form 

B k = [ dvol (z) try 1 3? □ k ' 2 ^t + J2 ® V ^ V 2fe - 6_i 9? 

Ij I 0<i<2fc-6 

+ •••+ ]T ^(V^)^~ i_3 (V5R) +U k \. (106) 

0<i<fc-3 J 



4.2 Leading derivatives in heat kernel asymptotics 

More precisely, all quadratic terms can be reduced to a finite number of invariant struc- 
tures, viz. m u 



+fi 3 ] Q □ k - 2 R + ff ] R^ □ k - 2 W u + fl 5) R □ fe - 2 i?| , (107) 

where R^ v is the Ricci tensor and fjp are some numerical coefficients. These numerical 
coefficients can be computed by the technique developed in the previous section. From 
the formula we have for the diagonal coefficients [bk] up to cubic terms in the jets 

[b k ] = <0\b k >= 1 <0;fc-l|L|0> 

k ) 

k-l 2(fc-i-l) pA-A 
+ (-l) fc E E TiCT (WL-i) <0;k-i-l\L\n i >< n t ; i - 1|L|0 > 



1=1 rij=0 I 



+0(K 3 ), (10* 



where 



< n;k\L\m >= (W k g 1 )-k<n\L\m> (109) 

and 0(3? 3 ) denote terms of third order in the jets. 

By computing the matrix elements in the second order in the jets and integrating over 
M one obtains 0], § 

= 1 (110) 

/( 2 ) 1 
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f ^ " 2(4P-1) (113) 

„r^ k 2 — k — 1 , 

ft = • (114) 

Jk 4(4P-1) 1 J 

One should note that the same results were obtained by a completely different method 
by Branson, Gilkey and 0rsted W2\. 



4.3 'Summation' of asymptotic expansion 

Let us consider the situation when the curvatures are small but rapidly varying, i.e. the 
derivatives of the curvatures are more important than the powers of them. Then the 
leading derivative terms in the heat kernel are the largest ones. Thus the trace of the 
heat kernel has the form 

Tr £,2 exp(-tF) ~ {47it)~ d/2 jfio - tB 1 + ^H 2 (t)\ + 0(K 3 ), (115) 

where H2(t) is some complicated nonlocal functional that has the following asymptotic 
expansion as t — » 

#a(t) ~ 2 E {j ^n— B ^ + °(^ 3 )- ( 116 ) 

k>2 

Using the results for Bk,2 one can easily construct such a functional H 2 just by a formal 
summing the leading derivatives 

H 2 = J dvol{x)tT V lQ^ 1) {-tn)Q + 2n afl V a ^ {2) {-tD)V u n^ 

-2Q 7 (3) (-t d)R + R^ A \-tn)R» u + R^X-t n)i?j, (117) 
where y®(z) are entire functions defined by |J1 0] 



i 



7W( " ) = ^ (2lTlj! / ^ = /^/ (i) (Oexp(-^^) (H8) 



o 



where 



/ (1) (0 = 


1 


/ (2) (0 = 


¥ 


/ (3) (0 = 




/ (4) (0 = 




/ {5) (0 = 


4l< 3 



-2 tf 4\ 



(119) 
(120) 

(121) 

(122) 

(123) 
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Therefore, H2(t) can be regarded as generating functional for quadratic terms B^p (leading 
derivative terms) in all HMDS-coefficients B^. It plays also a very important role in 
investigating the nonlocal structure of the effective action in quantum field theory in 
high-energy approximation 0. 

Let us note also that the function B q introduced in the Sect. 1 can be obtained just 
by analytical continuation of the formula for B^ from integer points A; to a complex plane 



4.4 Covariantly constant background 

Let us consider now the opposite case, when the curvatures are strong but slowly varying, 
i.e. the powers of the curvatures are more important than the derivatives of them. The 
main terms in this approximation are the terms without any covariant derivatives of the 
curvatures, i.e. the lowest order jets. We will consider mostly the zeroth order of this 
approximation which corresponds simply to covariantly constant background curvatures 

VRiem = 0, VU = 0, VQ = 0. (124) 

The asymptotic expansion of the trace of the heat kernel 

Tr L2 exp(-tF) ~ (47ri)- d / 2 ]T ^B k . (125) 

fe>0 fc - 

determines then all the terms without covariant derivatives (highest order terms in the 
jets), B^k, in all HMDS-coefficients B^. These terms do not contain any covariant deriva- 
tives and are just polynomials in the curvatures and the endomorphism Q. Thus the trace 
of the heat kernel is a generating functional for all HMDS-coefficients for a covariantly 
constant background, in particular, for all symmetric spaces. Thus the problem is to 
calculate the trace of the heat kernel for covariantly constant background. 



4.4.1 Algebraic approach 

There exist a very elegant indirect way to construct the heat kernel without solving the 
heat equation but using only the commutation relations of some covariant first order 
differential operators 0, [|, [5|, ||. The main idea is in a generalization of the usual Fourier 
transform to the case of operators and consists in the following. 

Let us consider for a moment a trivial case, where the curvatures vanish but not the 
potential term: 

Riem = 0, K = 0, VQ = 0. (126) 

In this case the operators of covariant derivatives obviously commute and form together 
with the potential term an Abelian Lie algebra 

[V„,V„] = 0, [V„,Q] = 0. (127) 

It is easy to show that the heat semigroup operator can be presented in the form 

exp(-tF) = (Ant)- d/2 exp(-tQ) J dvol (Jfe) exp (-— < k, gk > +k ■ v) , (128) 
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where < k, gk >= k^g^k", k- V = k^V^. Here, of course, it is assumed that the covariant 
derivatives commute also with the metric 

[V,g} = 0. (129) 

Acting with this operator on the Dirac distribution and using the obvious relation 

ex P (k-\7)6(x,x')\ x=x ,=6(k), (130) 

one integrates easily over k and obtains the trace of the heat kernel 

Tr L2 exp(-tF) = (47rt)~ d/2 J dvol (x)tr v exp(-tQ). (131) 

M 

In fact, the covariant differential operators V do not commute, their commutators 
being proportional to the curvatures 9ft. The commutators of covariant derivatives V 
with the curvatures 9ft give the first derivatives of the curvatures, i.e. the jets 3ft(i), the 
commutators of covariant derivatives with 9ft(i) give the second jets 3ft(2), etc. Thus the 
operators V together with the whole set of the jets J form an infinite dimensional Lie 
algebra £ = {V,%;(z = 1,2,...)}. 

To evaluate the heat kernel in the considered (low-energy) approximation one can take 
into account a finite number of low-order jets, i.e. the low-order covariant derivatives of the 
background fields, {9ft(i); (i < N)}, and neglect all the higher order jets, i.e. the covariant 
derivatives of higher orders, i.e. put 9ft(j) = for i > N. Then one can show that there 
exist a set of covariant differential operators that together with the background fields and 
their low-order derivatives generate a finite dimensional Lie algebra Gn = {V, (i = 
1,2,..., AT)} iii- 



Thus one can try to generalize the above idea in such a way that ( |128| ) would be 
the zeroth approximation in the commutators of the covariant derivatives, i.e. in the 
curvatures. Roughly speaking, we are going to find a representation of the heat semigroup 
operator in the form 

exp(-tF) = J dk$(t,k)exp(-— < k,V(t)k > +k-T^J (132) 

where < k, #(t)Jfe >= k A ^ AB (t)k B , k-T = k A T Al (A = 1, 2, . . . , D), T A = A^V M + Y A are 
some first order differential operators and the functions \&(t) and $(£, k) are expressed in 
terms of commutators of these operators — i.e., in terms of the curvatures. 

In general, the operators T A do not form a closed finite dimensional algebra because 
at each step taking more commutators there appear more and more derivatives of the 
curvatures. It is the low-energy reduction Q — ► Q N , i.e. the restriction to the low-order 
jets, that actually closes the algebra Q of the operators T A and the background jets, i.e. 
makes it finite dimensional. 

Using this representation one can, as above, act with exp(/c • T) on the Dirac dis- 
tribution to get the heat kernel. The main point of this idea is that it is much easier 
to calculate the action of the exponential of the first order operator k ■ T on the Dirac 
distribution than that of the exponential of the second order operator □. 
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4.4.2 Covariantly constant bundle curvature and covariantly constant endo- 
morphsim Q in flat space 

Let us consider now the more complicated case of nontrivial covariantly constant curvature 
of the connection on the vector bundle V in flat space: 

Riem = 0, VK = 0, VQ = 0. (133) 

Using the condition of covariant constancy of the curvatures one can show that in this 
case the covariant derivatives form a nilpotent Lie algebra H 

[V M , V„] = 7^, (134) 
[V A4J ^] = [V M ,Q] = J (135) 
[n^,n af3 ] = [TZ^,Q]=0. (136) 

For this algebra one can prove a theorem expressing the heat semigroup operator in 
terms of an average over the corresponding Lie group [[| 

exp(-tF) = (47rt)-^ exp(-tQ)det ^ 2 d(TM) ( ~ ^ ) (137) 
x J dvol (k) exp (-— < k, gtKcoth (ffi.)k > +k ■ , (138) 



d 



JR 

where k ■ V = PV^. Here functions of the curvatures 7£ are understood as functions 
of sections of the bundle End (TM) (g> End (V), and the determinant det End (tm) is taken 
with respect to End (TM) indices, End (V) indices being intact. 
It is not difficult to show that also in this case we have 

eMk-V)6(x,x')\ x=xl = 6(k). (139) 

Subsequently, the integral over fc M becomes trivial and we obtain immediately the trace 
of the heat kernel 

Tr i2 exp(-tF) = (4irt)- d / 2 J dvo\ (x)tr v exp(-tQ)det ^ d(TM) ^— ^—^ . (140) 

Expanding it in a power series in t one can find all covariantly constant terms in all 
HMDS-coefficients B k . 

As we have seen the contribution of the bundle curvature TZ^ is not as trivial as that 
of the potential term. However, the algebraic approach does work in this case too. It is a 
good example how one can get the heat kernel without solving any differential equations 
but using only the algebraic properties of the covariant derivatives. 



Ivan G. Avramidi: Computation of the heat kernel 



24 



4.4.3 Contribution of two first derivatives of the endomorphism Q 

In fact, in flat space it is possible to do a bit more, i.e. to calculate the contribution of the 
first and the second derivatives of the potential term Q [|J. That is we consider the case 
when the derivatives of the endomorphism Q vanish only starting from the third order, 
i.e. 

Riem = 0, VIZ = 0, VVVQ = 0. (141) 

Besides we assume the background to be Abelian, i.e. all the nonvanishing background 
quantities, TZ a /3, Q, Q-^ = V M Q and Q ]Ufl = V^V^Q, commute with each other. Thus we 
have again a nilpotent Lie algebra 

[V M , V,] = (142) 
[V M ,Q] = Q iAt (143) 
[V„, Q ; „] = (144) 

all other commutators being zero. 

Now, let us represent the endomorphism Q in the form 

Q = Q- a ik NiN k , (145) 

where (z = 1, . . . , q; q < d), a lk is some constant symmetric nondegenerate q x q matrix, Q 
is a covariantly constant endomorphism and iVj are some endomorphisms with vanishing 
second covariant derivative: 

Vfi = 0, VViV; = 0. (146) 
Next, let us introduce the operators Xa = (V M , iVj), (A = 1, . . . , d + q) and the matrix 

™ = (X. N T)' < 147 > 

with = V^. 

The operator F can now be written in the form 

F = -G AB X A X B + Q, (148) 

where 



(G")=^0 a *J- ("9) 
and the commutation relations (|144 ) take a more compact form 

\X a ,X b \=Tab (150) 



all other commutators being zero. 
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This algebra is again a nilpotent Lie algebra. Thus one can apply the previous theorem 
in this case too to get [§J 



exp(-tF) = (47rf) 



-(d+g)/2 



exp(— if2)det 



1/2 



tT 



x 



y sinh (t J 7 ) J 

J dkG 1/2 exp (-^ < k, GtFcoth {tT)k > +k ■ x) , (151) 



JR d +i 



where G = det Gab and k ■ X — k A X&- 

Thus we have expressed the heat semigroup operator in terms of the operator exp(k ■ X) . 
The integration over k is Gaussian except for the noncommutative part. Splitting the in- 
tegration variables (k A ) = (q^,u z ) and using the Campbell-Hausdorf formula we obtain 



exp(A; • X)5(x, x ) _ / = exp(co> • N)5(q), 



(152) 



where uo ■ N = uj l Ni. 

Further, after taking off the trivial integration over q and a Gaussian integral over u, 
we obtain the trace of the heat kernel in a very simple form M 



Tr L 2 exp(-tF) = (4vrf)~ d/2 J dvo\ (x)tr v $(t) exp 



M 



-tQ + -t 3 <VQ,^(t)g- 1 VQ> 



where < VQ, V^g^VQ >= V ^Q^ u {t)g uX V X Q 

\ + t 2 C{t)P 



(153) 



det ^% M) K{t) det ~H d \ TM) {l + t 2 [E(t) - S(t)K-\t)S(t)]P} 



x det" 1/2 

AUt,L End (TM) 



*(*) = {*"„(*)} = [i+t 2 c(t)py c(t), 

P is the matrix determined by second derivatives of the potential term, 

P = {P*} , Pt = I^v,v a q, 



(154) 
(155) 

(156) 



and the matrices C{t) = {C^{t)}, K(t) = {K» u {t)} S(t) = {S%{t)} and E(t) = {E^{t)} 
are defined by 

C(t)= l — tcot\i(tz- l )(l- zK- z 2 P)-\ (157) 
J 2ni 

c 



K(t) 



c 



dz t 
2m z 2 



S(t) = I - sinh -zTZ- z 2 P)-\ 

J 2m z 

c 



sinh (tz~ x ){l -zK- z 2 P)-\ 



(158) 
(159) 
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E(t)= I ^ : tsmh(^- 1 )(l - zK- z 2 P)-\ (160) 
J 2m 

c 

where the integral is taken along a sufficiently small closed contour C that encircles the 
origin counter-clockwise, so that F(z) = (1 — zTZ — z 2 P)~ 1 is analytic inside this contour. 

The formula ( |153| ) exhibits the general structure of the trace of the heat kernel. 
Namely, one sees immediately how the endomorphism Q and its first derivatives VQ 
enter the result. The nontrivial information is contained only in a scalar, and a 

tensor, ty^ u (t), functions which are constructed purely from the curvature TZ^ U and the 
second derivatives of the endomorphism Q, VVQ. 

So, we conclude that the HMDS-coefficients Bk are constructed from three different 
types of scalar (connected) blocks, Q, $( n) (ft, VVQ) and V^Qtfj^ (ft, VVQ)V„Q. They 
are listed explicitly up to B 8 in N7J. 



4.4.4 Symmetric spaces 

Let us now generalize the algebraic approach to the case of the curved manifolds with 
covariantly constant Riemann curvature and the trivial bundle connection H || 

VRiem = 0, TZ = 0, VQ = 0. (161) 

First of all, we give some definitions. The condition ( |161| ) defines, as we already 
said above, the geometry of locally symmetric spaces. A Riemannian locally symmetric 
space which is simply connected and complete is globally symmetric space (or, simply, 
symmetric space) . A symmetric space is said to be of compact, noncompact or Euclidean 
type if all sectional curvatures K(u, v) = R a b c dU a v b u c v d are positive, negative or zero. A 
direct product of symmetric spaces of compact and noncompact types is called semisimple 
symmetric space. A generic complete simply connected Riemannian symmetric space is 
a direct product of a flat space and a semisimple symmetric space. 

It should be noted that our analysis in this paper is purely local. We are looking 
for a universal local function of the curvature invariants, that reproduces adequately the 
asymptotic expansion of the trace of the heat kernel. This function should give all the 
terms without covariant derivatives of the curvature in the asymptotic expansion of the 
heat kernel, i.e. in other words all HMDS-coefficients Bk for any locally symmetric space. 

It is well known that the HMDS-coefficients have a universal structure, i.e. they are 
polynomials in the background jets (just in curvatures in case of symmetric spaces) with 
the numerical coefficients that do not depend on the global properties of the manifold, 
on the dimension, on the signature of the metric etc. It is this universal structure we are 
going to study. 

It is obvious that any flat subspaces do not contribute to the HMDS-coefficients Bk- 
Therefore, to find this universal structure it is sufficient to consider only semisimple 
symmetric spaces. Moreover, since HMDS-coefficients are analytic in the curvatures, one 
can restrict oneself only to symmetric spaces of compact type. Using the factorization 
property of the heat kernel and the duality between compact and noncompact symmetric 
spaces one can obtain then the results for the general case by analytical continuation. 
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That is why we consider only the case of compact symmetric spaces when the sectional 
curvatures and the metric are positive definite. 

Let e a be a basis in the tangent bundle which is covariantly constant (parallel) along 
the geodesic. The frame components of the curvature tensor of a symmetric space are, 
obviously, constant and can be presented in the form 

Rated = PikE\ b E k cd , (162) 

where E l ab , (i = 1, . . . ,p;p < d(d— l)/2), is some set of antisymmetric matrices and f3i k is 
some symmetric nondegenerate p x p matrix. The traceless matrices Di = {D a ib } defined 
by 

D\ b = -ftfc^V = ~ D \i ( 163 ) 
are known to be the generators of the holonomy algebra H, 

[Di, D k ] — F j ik Dj, (164) 

where F\ k are the structure constants. 

In symmetric spaces a much richer algebraic structure exists. Indeed, let us define the 
quantities C A BC = -C A OB , (A = 1, . . . , D; D = d + p): 

C\b = F l ab , C a ib = D a ib , C l kl = F l kl , (165) 

C\ c = C\ a = G\ = 0, (166) 
and the matrices Ca = {C B AC } = (C a ,Ci): 

C=(£ C,^(^ »V (1 67) 



\EL J' 10 F ik 

One can show that they satisfy the Jacobi identities 0, ^| 

[C A ,C B ] = C C AB C C (168) 

and, hence, define a Lie algebra Q of dimension D with the structure constants C A BC , the 
matrices C A being the generators of adjoint representation. 

In symmetric spaces one can find explicitly the generators of the infinitesimal isome- 
tries, i.e. the Killing vector fields £4, and show that they form a Lie algebra of isometries 
that is (in case of semisimple symmetric space) isomorphic to the Lie algebra Q, viz. 

[UM=C° AB ^c. (169) 
Moreover, introducing a symmetric nondegenerate D x D matrix 

7AB=f? iM, (170) 



V A 



ik 



that plays the role of the metric on the algebra Q, one can express the operator F in 
semisimple symmetric spaces in terms of the generators of isometries 



F = - 1 AB UiB + Q, (171) 
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where j AB = (jab)" 1 - 

Using this representation one can prove a theorem that presents the heat semigroup 
operator in terms of some average over the group of isometries G H || 



exp(-iF) = (47rf)~ D/2 exp 



t[Q- -R G 



(172) 



x 



dkj 1/2 det l H {g) 



/ sinh(fc-C/2) \ 
^ k-C/2 ) 



exp 



-—<k^k>+k-i 



where 7 = det7AB, k ■ C = k A CA, fc • £ = k A ^A, and Rg is the scalar curvature of the 
group of isometries G 



R 



G 



^7 ^ AD° BC- 



(173) 



Acting with this operator on the Dirac distribution 5(x,x f ) one can, in principle, 
evaluate the off-diagonal heat kernel exp(— tF)8(x, x'), i.e. for non-coinciding points x 7^ 
x' (see 0). To calculate the trace of the heat kernel, it is sufficient to compute only the 
coincidence limit x = x'. Splitting the integration variables k A = (q a ) uo' 1 ) and solving 
the equations of characteristics one can obtain the action of the isometries on the Dirac 
distribution || |6| 



exp (k ■ £) 5(x, x') 



utiL End (TM) 



/sinh {u-D/2y 



S(q). 



(174) 



where uj ■ D = uj l Di. 

Using this result one can easily integrate over q in (|172 ) to get the heat kernel diagonal. 
After changing the integration variables uj — > y/tuj it takes the form || 



P(t)) 



(Airt) 



-d/2 



exp 



—t ( Q R Ru 

8 6 



x(47r)- p / 2 J du(3 1/2 exp (-- <uj,puj> 



Ad(n) V y/iu ■ F/2 J 



det " 1/2 

ueL End (TM) 



fsmh(Viuj ■ D/2)\ 
V x/tuj-D/2 J 



(175) 



where uj ■ F = uJ l Fi, -Fj = {-F^} are the generators of the holonomy algebra Ti in adjoint 
representation and 

R H = ~(3 ik F^F l km (176) 



is the scalar curvature of the holonomy group. 

The remaining integration over uj in (|175 ) can be done in a rather formal way || |9[] . 
Let a* and a k be operators acting on a Hilbert space, that form the following Lie algebra 



\a J ,a 



5 "fcj 



k ■ 



(177) 
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[a\a k ] = [a*,a* k ] = 0. 



(178) 



Let |0 > be the 'vacuum vector' in the Hilbert space, i.e. 



< 0|0 >= 1, 



(179) 



a*|0 >= 0, < 0\a* k = 0. (180) 

Then the heat kernel ( |175| ) can be presented in an formal algebraic form without any 
integration 





xexp(< a*,/TV >) |o). (181) 



where a-F = a k F k and a-D = a k D k . This formal solution should be understood as a power 
series in the operators a k and a* k and determines a well defined asymptotic expansion in 



Let us stress that these formulas are manifestly covariant because they are expressed 
in terms of the invariants of the holonomy group H, i.e. the invariants of the Riemann 
curvature tensor. They can be used now to generate all HMDS-coefflcients [b k ] for any 
locally symmetric space, i.e. for any manifold with covariantly constant curvature, simply 
by expanding it in an asymptotic power series as t — > 0. Thereby one finds all covariantly 
constant terms in all HMDS-coefflcients in a manifestly covariant way. This gives a very 
nontrivial example how the heat kernel can be constructed using only the Lie algebra of 
isometries of the symmetric space. 

5 Conclusion 

In present paper we have presented recent results in studying the heat kernel obtained 
in our papers [Jj, ^, ^ f|, |5|, |5[]. We discussed some ideas connected with the problem of 
developing consistent covariant approximation schemes for calculating the heat kernel. 
Especial attention is payed to the low-energy approxiamtion. It is shown that in the local 
analysis there exists an algebraic structure (the Lie algebra of background jets) that turns 
out to be extremely useful for the study of the low-energy approximation. Based on the 
background jets algebra we have proposed a new promising approach for calculating the 
low-energy heat kernel. 

Within this framework we have obtained closed formulas for the heat kernel diagonal 
in the case of covariantly constant background. Besides, we were able to take into ac- 
count the first and second derivatives of the endomorphism Q in flat space. The obtained 
formulas are manifestly covariant and applicable for a generic covariantly constant back- 
ground. This enables to treat the results as the generating functions for the whole set of 



t 



0. 
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the Hadamard-Minakshisundaram-De Witt-Seeley-coefficients. In other words, we have 
calculated all covariantly constant terms in all HMDS-coefficients. 

Needless to say that the investigation of the low-energy effective action is of great 
importance in quantum gravity and gauge theories because it describes the dynamics of 
the vacuum state of the theory 
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